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Abstract. - The quantum dynamics of two entangled two-level atoms is studied. Each of the two 
atoms is located within an isolated and dissipative cavity. If the interaction time of atoms and 
cavities is not very long, the amount of quantum discord and entanglement between two atoms 
decreases as the system evolves. The sudden death of quantum discord and entanglement of two 
atoms occurs within a short interaction time. However, after a long interaction time, quantum 
discord and entanglement of two atoms could be partially preserved due to the long-lived nature 
of quantum discord and entanglement. Surprisingly, we find the amount of long-lived quantum 
discord could be smaller than that of long-lived entanglement. Thus, entanglement may be more 
robust than quantum discord against decoherence. 



Introduction. — Quantum entanglement is at the 
heart of quantum information processing and quantum 
computation [TH3]. In recent years, many efforts have 
been invested in the study of the evolution of joint sys- 
tems formed by two subsystems (each system locally in- 
teracts with its environment) [IHS]- In particular, the en- 
tanglement of a two-qubit system may disappear for a fi- 
nite time during the dynamics evolution. The nonsmooth 
finite-time disappearance of entanglement is called "entan- 
glement sudden death" (ESD). Experimentally, the BSD 
phenomenon has been observed in the laboratory by sev- 
eral groups for optical setups p^lfTT] and atomic ensembles 

One the other hand, quantum entanglement is not the 
only kind of quantum correlation useful for quantum in- 
formation processing [T3HT5] . In fact, it was shown both 
theoretically [I6l - l23] and experimentally [21] that some 
tasks can be sped up over their classical counterparts us- 
ing fully separable and highly mixed states. These results 
clearly show that separable states with quantum discord 
can be used to implement quantum information process- 
ing such as deterministic quantum computation with one 



qubit (24j . Quantum discord introduced in |25 1 I26 | is an- 
other kind of quantum correlation different from entan- 
glement. Very recently, quantum discord has been inves- 
tigated widely [27^150] . Note that all the previous stud- 
ies p7H30] have shown that, for several quantum systems, 
there is no quantum discord sudden death (DSD). How- 
ever, quantum discord is a kind of quantum correlation in 
composite quantum systems. Since entanglement of quan- 
tum systems can stay zero for a finite time (ESD), a natu- 
ral question is whether there is DSD in quantum systems 
with ESD. Here, we present a quantum system where there 
exists DSD as well as ESD. We also explain why there is 
no DSD in [3^. Furthermore, we find that there is also 
long-lived quantum discord. In recent years, many efforts 
has been devoted to the study of the long-lived entangle- 
ment in cavity QED [3TH34] or solid state systems [35] . To 
the best of our knowledge, there is few study on the long- 
term behavior of quantum discord. Thus, an investigation 
of quantum discord of a quantum system in the presence 
of decoherence in the limit i — >■ cxi is highly desired. This 
question is also addressed in the present work. 

In the present paper, we investigate the dynamics of 
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quantum discord and entanglement of a quantum system 
formed by two two-level atoms within two spatially sepa- 
rated and dissipative cavities in the dispersive limit using 
the results of |36| . The two atoms are initially prepared 
in the Werner states [5] and the cavities are initially pre- 
pared in coherent states. We show that both DSD and 
ESD can appear in the present system. The amount of 
quantum discord and entanglement of two atoms decreases 
with time in the short-term. However, the long-term be- 
havior is very different since a long survival of quantum 
discord and entanglement are shown in the system. This 
implies that quantum discord and entanglement of two 
atoms could be partially preserved even they are put into 
dissipative cavities. Unlike the results in [27l|30], our re- 
sults show that the amount of long-lived quantum discord 
could be smaller than that of long-lived entanglement. In 
other words, quantum entanglement may be more robust 
than quantum discord in the present model. 

The model. — We first consider a quantum system 
consisting of a two-level atom interacting with a single- 
mode cavity. Under the electric dipole and rotating wave 
approximation, the Hamiltonian of the present system is 

{h = i)m 



H 



(1) 



where g is the atom- field coupling constant, a± are the 
atomic spin flip operators characterizing the effective two- 
level atom with frequency wqi and cr^ — |e)(e| — |5)(<?|- 
Note that the symbols |e) and \g) refer to the excited and 
ground states for the two-level atom. Here, a) and a are 
the creation and annihilation operators of the field with 
frequency uj, respectively. The dispersive limit is obtained 
when the condition |A| = jwo — wj 3> \/n+ Ig is satisfied 
for any relevant n. Then, the interaction Hamiltonian 
g{a^a- + aa+) can be regarded as a small perturbation. 
Hence the effective Hamiltonian of the present model can 
be rewritten as j38) 



iJe = wa^a + -^a^ + n[{a^a + l)|e)(e| - a^a\g){g\], 



(2) 



with Q = (?^/A. In the interaction picture, the interaction 
Hamiltonian is 



V ^n[{a^a + l)\e){e\-a^a\g){g\]. 



(3) 



We assume the two-level atom interacting with a co- 
herent field in a dissipative environment. This interaction 
causes the losses in the cavity which is presented by the 
superoperator D = 7 (2a ■ — a^a ■ — ■ a^a), where 7 is 
the decay constant. For the sake of simplicity, we confine 
our consideration in the case of zero temperature cavity. 
Then, the master equation that governs the dynamics of 
the system can be written as follows 



where p is the density matrix of the atom-field system. 

Ca Cc\ 

C Cb J 

and the field is initially prepared in a coherent state 



If the initial state of the two-level atom is 



=0 v^i 



with a being a complex num- 
n\n). Then, 



ber. Here, \n) is the Fock state with a' a|n) = rt 
the reduced density matrix of the atom is obtained by 
tracing out the variables of the field from the atom-field 
density matrix 1391 



Patom (^) 
fit) 



Ca\e){e\+(:t\g){g\ + [Ccf{t)\e){g\ + h.c], 
exp{-im+ |a|2(e-27t _ 1)} 

-2(7+jO)tj| 

1)}, (5) 



X exp{ 1 
xexp{|ape-2''*(e 



where h.c denotes the Hermitian conjugate. 

Then, we consider a quantum system consisting of two 
noninteracting atoms each locally interacts with its own 
coherent field of a dissipative cavity. The interactions be- 
tween each atom and its own dissipative cavity is described 
by Eq. (|4]). We assume the two atoms are initially pre- 
pared in Werner states defined by '6' 

p$ = p|ci>)(ci>| + i_P/, 

p^ = p|vI/)(vI/| + i_Z/, 

I*) = -^i\eg) + \ge)), 

l*> = -^i\ee) + \gg)), (6) 

where p is a real number which indicates the purity of 
initial states, / is a 4 x 4 identity matrix. The parameter 
p is I for pure sates and for completely mixed states. 
The two fields are prepared in coherent states |ai) and 
|q!2)- For the sake of simplicity, we assume ai = a2 = a, 
the decay rates of the two cavities are equal, and the atom- 
field coupling constants are the same. Using the method 
introduced in , we can obtain the reduced density matrix 
of two atoms conveniently. 

The reduced density matrix of two atoms can be ob- 
tained by using the superoperator method [39ti4Tj . As one 
can see below, the quantum discord and entanglement of 
the present system can be calculated conveniently by em- 
ploying the results of [3^. We assume the initial state of 
two atoms is p^. Using the results of [39] and Eq. ([6]), we 
obtain the density matrix of two atoms p(t) as follow 



V 
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(7) 



(4) where f{t) is given by Eq. ([5]). 
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Quantum discord and entanglement. — In gen- 
eral, a composite quantum system contains both quantum 
and classical correlations, the total amount of which are 
quantified by quantum mutual information. Precisely, the 
quantum mutual information of a composite bipartite sys- 
tem p^^ is defined as 



AB\ 



Sip"") + S{p") - ^(p^^). 



(8) 



where S{p) = —Tr{p\og2 p) = — lo§2 ^i) '^^n 

Neumann entropy of density matrix p with being the 
eigenvalues of density matrix p. We note that logj is 
defined to be and p^{p^) is the reduced density matrix 
of p"^^ by tracing out system B[A). 

Quantum discord 25, 26^ is another kind of quantum 
correlation different from entanglement. In order to quan- 
tify quantum discord, the authors of [25* proposed to use 
the von Neumann type measurements consisting of one- 
dimensional projector {Bi\ (acts on system B only), such 
that "^^Bi = 1. The conditional density matrix of the 
total system after the von Neumann type measurements 
is [Ml 



1 



{I®Bi)p''''{I®B.,), 



p, = Tr{{I®Bi)p'^"{I®Bi)) 



(9) 



where Pi is the probability of the corresponding measure- 
ment. The quantum conditional entropy with respect to 
this kind of measurement is defined as 



s{p^^m})^Y.p^s{pii 



(10) 



and the corresponding quantum mutual information with 
respect to the measurement is defined by 



Z(p^^|{S,})^5(p-^)-^(p-«|{6.}) 



AB\ 



(11) 



The quantity I{p'^^\{Bi\) is the information gained about 
system A if one performs measurement Bi on system B. 
The resulting classical correlation according to [25jl26j is 
defined as 

j{p^^) = supz(p^^|{sa) 

= 5(/)-min[5(/^|{S.})]. (12) 

The quantum discord is obtained by subtracting J from 
the quantum mutual information X 



2?(p^^)=Z(p^^)-^(p^^) 



(13) 



As one can see from the above equations, the minimization 
procedure should be done over all possible von Neumann 
measurements Bi on system B. Thus, the main difficulty 
of calculating quantum discord lies in the elaborate mini- 
mization process in the term miniej [<S'(/3'^^|{Si})]- Fo^" 
tunately, the quantum discord of Eq. ([7]) can be calculated 



with the help of the results of 
1 



V{p 



AB\ 



-[(l-rfi-d2 

+ {l-di +d2 

-|-(l + fii -d2 
-|-(l + rfl +d2 
1 



d , 1-d 
-log2 



2 °^ 2 2 
di = d2^p\f{t)\^,d3 = -p, 

d = max{|di|, 1^21, Mall- 



- c?3)log2(l - di - 

c?3)log2(l -rfi 4 
c?3)log2(l + di- 

fi3)l0g2(l + rfl 4 

1 + rf, 1 

l0g2 - 



d2 
d2 
d2 



d3)] 
d3)] 
d3)] 



(14) 



In order to investigate the entanglement of two-qubit 
systems, we adopt the entanglement measure concurrence 
introduced in 1421 



C = max{0,xi - X2 - X3 - Xi}, 



(15) 



where Xi (* — 17 2,3,4) are the square roots of the eigen- 
values in decreasing order of the magnitude of the "spin- 
flipped" density matrix operator i? — p{ay®(Ty)p* {(Ty®ay) 

and tTj, is the Pauli Y matrix, i.e., Oy — (^^^ O* ) ' 

currence of a quantum state ranges from 0, which corre- 
sponds to an unentangied state, to 1, which corresponds 
to a maximally entangled state. The concurrence of the 
above state is 



1 



C(t)=max{0,p|/(i)|'-^} 



(16) 



We will use this equation to calculate the entanglement of 
the quantum system presented in this work. 

Results and discussions. — 

Sudden death of quantum discord and entanglement. 
We now want to investigate the dynamics of quantum dis- 
cord and entanglement. In Fig. [TJ quantum discord and 
entanglement are plotted as functions of the dimensionless 
scaled time Q,t for a = 0.5 (upper panel) and a = 1 (lower 
panel). Clearly, there is ESD in the present model as one 
can see from Fig. [TJ From the upper panel of Fig. [l] one 
may conclude that there is no DSD, which is consistent 
with the results of [271 [30]. However, this is not always 
correct. As one can easily observe from the lower panel of 
Fig. [U if we increase the intensity of the coherent fields 
that is proportional to jap, there is DSD. We note that in 
this figure, quantum discord is larger than entanglement, 
which is coincidence with the observations of \27_30 . How- 
ever, this is not a general result since quantum discord and 
entanglement are two different quantities and there is no 
simple relative ordering between them. For example, for 
the Werner states defined by p = p|V'-)(V'-| + where 
I ■(/'_) = {\eg) — \ge))/^/2, quantum discord may be lager 
or smaller than entanglement [36j . 

We now want to explain why there is no DSD in the case 
of dephasing channel in 00]. At first sight, the elements 
of the density matrix of the dephasing case is very similar 
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Fig. 1: Quantum discord (solid line) and entanglement (dotted 
line) of two atoms are plotted as functions of the dimensionless 
scaled time Qt with 7/f2 — 0.01 and p — 0.4. Upper panel: 
a = 0.5. Lower panel: a — 1. There is no DSD in the upper 
panel with a — 0.5. However, the lower panel clearly shows 
the existence of DSD in the case of a = 1. 

to the matrix elements of Eq. ([7]). Let us focus on the off- 
diagonal elements P23(t) of [20], i-e., P23{t) = e^^*p23(0), 
where T is the decay rate. Obviously, the term e~^* be- 
comes zero only in the asymptotic limit t — >■ oo. Quantum 
discord vanishes only in the asymptotic limit, which be- 
haves similarly to decoherence of each atom [30]. In the 
present work, the off-diagonal element P23(t) is much more 
complicated than that of 30 and it is possible for quan- 
tum discord to stay zero for a finite time. Physically, the 
infiuence of the interactions between atoms and cavities, 
and the mean photon number and decay rate of dissipa- 
tive cavities upon quantum discord has not bee considered 
in [30) . Here, all the above influence upon the dynamics 
of quantum discord and entanglement is taken into ac- 
counted. In this sense, the results of [301 can only re- 
veal parts of the properties of quantum discord (no DSD), 
which is consistent with the upper panel of Fig. [1] How- 
ever, the lower panel of Fig. [TJ which indicates the exis- 
tence of DSD, can not appear in [30 . 

Long-lived quantum discord and entanglement. As we 
have pointed out previously, little attention has been paid 
to long-lived quantum discord even though lots of work 
has been made on long-lived entanglement [3TH35] . Here, 
we consider the question of whether there is long-lived 
quantum discord and entanglement in the present model. 
The influence of the purity of the initial state of atoms 
upon the long-time behavior of quantum discord and en- 
tanglement is also discussed. 

In Fig. [21 quantum discord and entanglement as func- 
tions of the dimensionless scaled time ^t are plotted for 
p = 0.5 (upper panel) and p = 0.8 (lower panel). It is 
easy to observe that the long-time behavior of quantum 
discord and entanglement depends on the purity of the ini- 
tial state of atoms. This figure is a direct evidence of the 
presence of long-lived quantum discord and entanglement 
in the present system. Comparing the upper panel with 
the lower one of Fig. [2] one can see that the amount of 
long-lived quantum discord and entanglement will increase 
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Fig. 2: Quantum discord (solid line) and entanglement (dotted 
line) of two atoms are plotted as functions of the dimensionless 
scaled time Qt with 7/(1 = 0.01, and a = 0.5. Upper panel: 
p = 0.5. Lower panel: p = 0.8. Comparing the upper and lower 
panel, we see that long-lived quantum discord could be larger 
(upper panel) or smaller (lower panel) than long-lived entan- 
glement. In other words, quantum entanglement may be more 
robust than quantum discord in the presence of decoherence. 
This result is different from the previous work |27»30) 

with the increase of the parameter p. The relative ordering 
of quantum discord and entanglement depends heavily on 
the purity p. For example, in the case oip = 0.5, long-lived 
quantum discord is larger than long-lived entanglement. 
However, in the case of p = 0.8, we find that, in contrast 
to the results of [?7ll5n] , the amount of long-lived quantum 
discord could be smaller than that of long-lived entangle- 
ment. In other words, quantum entanglement could be 
more robust than discord against decoherence. Thus, it 
is difficult for us to make a general statement whether or 
not quantum discord is more robust against decoherence 
than entanglement. Intuitively, quantum discord is differ- 
ent from entanglement and it is difficult to make a general 
conclusion about the relative ordering of the amount of 
quantum discord and entanglement in the presence of de- 
coherence. In order to show the influence of decoherence of 
cavities upon the quantum discord of two atoms, we plot 
the quantum discord of two atoms in Fig. 3. Note that it 
has been proved that white noise of cavity fields can play 
a constructive role in the generation of entanglement in 
cavity QED systems [J^. Here, Fig. [3] demonstrates that 
the dissipation of cavity fields may also play a construc- 
tive role in the generation of quantum discord. Note that 
in Figs. ([J) -([SI, we have assumed the atoms are initially 
prepared in p^. One can also consider the dynamics of 
quantum discord and entanglement if the initial state of 
atoms is and the results are similar. 

Conclusions. — In the present work, we have stud- 
ied the dynamics of quantum discord and entanglement 
of two two-level atoms without direct interactions. Each 
atom is put into a spatially separated and dissipative cav- 
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Fig. 3: Quantum discord of two atoms are plotted as functions 
of the dimensionless scaled time fit and parameter 7/^ with 
Q = 1 and p — 0.8. 

ity in the dispersive hmit. We first investigated the short- 
time behavior of quantum discord and entanglement in 
the presence of dissipation. The amount of quantum dis- 
cord and entanglement of two atoms decreases with time 
if the interaction time is not very long. Particularly, we 
have shown that both DSD and ESD could appear simul- 
taneously in the present system. This is different from the 
results of [371130]. Then, we discussed the long-time be- 
havior of quantum discord and entanglement of two atoms. 
We show there is long-lived quantum discord and entan- 
glement in the presence of the dissipation of cavities. Our 
results show that quantum entanglement may be more ro- 
bust against decoherence than quantum discord. Finally, 
we would like to point out that it could be possible to 
study quantum discord of any dimensional bipartite states 
from a geometrical point of view [33] . It is also interesting 
to compare the results of any dimensional bipartite states 
in the presence of decoherence with our work. 
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